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ABSTRACT 

The  Middleton  Clan  A  narrowband  non-Gaui*ian 
noiae  model  ia  conaidered.  The  performance  of  varioua 
simple  auboptimum  threshold  detectors  is  examined 
and  compared  with  the  performance  of  the  optimum 
detector.  For  example,  it  is  shown  that  for  most  cases 
of  interest  the  blanker  has  nearly  optimal  perfor¬ 
mance  while  the  soft  limiter  and  hard  limiter  have  sub¬ 
stantially  less  than  optimal  performance.  Then,  a  new 
approximation  to  the  locally  optimum  detector  non¬ 
linearity  is  developed  using  ideas  from  the  theory  of 
robust  detection. 

I.  INTRODUCTION 

Various  attempts  have  been  made  to  develop 
nodels  of  non-Gaussian  noises  (see,  for  example,  [1- 
3J).  Among  the  most  general  models  are  those 
leveloped  by  Middleton  [7-12],  Middleton  divides  non- 
•aussian  noise  into  two  classes.  A  and  B.  (There  has 
Also  been  consideration  of  a  Class  C  which  contains 
noises  which  are  sums  of  Class  A  and  Class  B  com¬ 
ponents  [9].)  Class  B  noises  are  broadband,  i.e.  those 
with  spectra  broader  than  the  passband  of  the 
receiver.  Class  A  noises  are  narrowband,  i.e.  have 
spectra  comparable  to  or  narrower  than  the  receiver 
passband. 

In  a  previous  paper  [l]  we  examined  the  Middleton 
Class  A  noise  model.  We  showed  that,  in  a  wide  variety 
of  cases,  the  first-order  noise  probability  density  func¬ 
tion  (PDF),  which  is  an  infinite  weighted  sum  of  Gaus¬ 
sian  PDF's,  can  be  closely  approximated  by  the  (nor¬ 
malized)  sum  of  the  first  two  terms  in  the  series.  We 
also  showed  that  the  detector  which  is  locally  optimum 
for  this  approximation  performs  very  well  for  the  origi¬ 
nal  Class  A  model. 

In  this  paper,  after  a  brief  review  of  previous  work 
in  Section  II.  we  examine  the  performance  of  other 
even  simpler  suboptimum  threshold  detectors.  We 
show  that  in  almost  every  case  at  least  one  of  these 
very  simple  detectors  is  nearly  optimal.  In  the  Section 
IV  we  return  to  the  two  term  approximation  discussed 
above.  Anticipating  that  this  simple  approximation 
would  not  work  in  every  case  we  developed  a  somewhat 
more  sophisticated  scheme  involving  ideas  from  robust 
detection.  This  scheme  is  discussed  in  Section  IV. 
n.  PREVIOUS  WORK 

In  order  to  formulate  his  model,  Middleton  [7-1 8] 
assumes  the  noise  has  the  form  X(t)  +  N(t)  where 
IV (f )  is  a  Gaussian  background  component  and 

xU) »  £  v,u.*)  (i) 
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where  Uf  denotes  the  received  waveform  from  an 
interference  source  and  d  is  a  random  parameter.  He 
then  assumes  that  the  waveform  receptions  are  Pois¬ 
son  distributed  in  time  and  shows  that  the  normalized 
(to  unit  variance)  noise  density  /  (x )  can  be  approxi¬ 
mated  canpnically  by 


/(*>*  t  Kmfl*:oZ)  (2) 
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where  f  (z  :e*)  is  the  zero-mean  Gaussian  PDF  with  vari¬ 
ance  «r.  The  variance  e£  of  the  m**  density  is  given  by 


i  _  m/A  +  P 
t  +  r 

and  the  coefficient  Km  is  given  by 


(3) 

(4) 


where  A  and  P  are  the  two  basic  parameters  of  the 
model.  The  first  parameter.  A.  is  called  the  "overlap 
index"  and  is  defined  by  A  =  vT.  where  v  is  the  rate  of 
the  homogeneous  Poisson  process  which  governs  the 
generation  of  the  interfering  waveforms  Ut  and  T  is  the 
mean  duration  of  a  typical  interfering  signal.  The 
other  parameter.  P,  is  given  by  the  ratio  or  the  power 
in  the  Gaussian  portion  of  the  interference  to  the 
power  in  the  Poisson  component. 

Middleton  has  shown  that,  by  adjusting  the  param¬ 
eters  A  and  P.  the  density  f  given  in  (2)  can  be  made  to 
fit  a  great  variety  of  non-Gaussian  noise  densities  [9- 
12].  Also,  the  parameters  A  and  P  are  physically 
motivated  and  can  be  directly  estimated  (see  [12,9]). 
Unfortunately  the  model  (2)  ia  cumbersome.  For 
example,  in  [11],  Spaulding  and  Middleton  exhibit  the 
optimal  nonlinearity  for  detection  (i.e.  the  likelihood 
ratio  /(*-*i)//(*-*o))  and  point  out  that  this  detec¬ 
tor  structure  is  likely  to  be  computationally  burden¬ 
some  and  uneconomical.  Thus  we  would  like  to  develop 
detector  nonlinearities  having  simpler  structure  but 
which  retain  the  desirable  properties  of  the  one  given 
in  [11]. 

In  [l]  we  conaidered  the  (normalized)  M-term 
truncation  of  the  Class  A  noise  PDF  given  in  (2).  i.e. 

t'  *./(*:<£) 

/(»>(.)  *  -  (5) 
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We  showed  that  /  <**  is  a  very  good  approximation  to  / 
tor  small  values  of  M.  In  particular,  for  the  problem  of 
threshold  (or  locally  optimum)  detection  (i.e.  small 
signal,  large  time-bandwidth  product),  the  detector 
designed  assuming  the  noise  has  density  /(<r>  performs 
virtually  as  well  as  the  optimal  detector  when  the  noise 
is  actually  given  by  the  full  Class  A  model  / .  even  for 
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In  the  next  section  we  will  examine  the  perfor¬ 
mance  of  several  suboptimum  detectors  whose  struc¬ 
ture  is  even  simpler  than  the  ones  considered  In  [l], 
but  first  we  will  give  a  brief  summary  of  the  necessary 
results  from  the  theory  of  locally  optimum  detection. 
The  details  of  this  theory  may  be  found  in  many  places 
[14-17.11], 


Under  mild  regularity  conditions  the  (asymptotic) 
performance  (or  processing  gain)  achievable  using  a 
given  detector  nonlinearity  g(z)  when  the  (indepen¬ 
dent  identically  distributed)  noise  process  has  first 
order  PDF  A(x)  is  given  by  the  efficacy  functional 
.  s 

f  g(z)K(z)dx 

n(g,h)  =  % - 1—  (6) 


/p*(x)A(x)dx 


For  a  given  noise  PDF  A(x)  the  locally  optimum  detec¬ 
tor  nonlinearity  gj(z )  is  the  solution  to 

ff  *  (M  =  max  r)(g,K)  (7) 

» 

and  has  the  form 

<» 

The  performance  of  this  locally  optimum  nonlinearity 
is  given  by 

»?•(*)  =  n(fv*)  =  f 

— m 

The  functional  tj*(A)  is  also  known  aa  Fisher's  Informa¬ 
tion  for  A.  An  interesting  (and  well  known)  fact  about 
the  function  t) *(A)  is  that  it  is  minimized  over  all  PDF's 
by  the  Gaussian  PDF.  In  fact,  the  locally  optimum 
detector  nonlinearity  for  Gaussian  noise  is  the  linear 
detector  (g  (*)=*/ or)  which  has  performance  equal  to 
unity  for  all  noise  PDF's. 

The  Class  A  noise  PDF  /  in  (2)  is  highly  non- 
Caussian,  and  it  is  often  the  case  that  t;*(/)>>1.  In 
[1]  it  was  shown  that  for  very  small  values  of  M,  /!*> 
and  gi(z)  closely  approximate  /  and  its  locally 
optimum  nonlinearity  gj(z).  Moreover,  for  a  wide 
range  of  values  of  the  parameters)  we  saw  that  the 
processing  gain  achievable  using  gi  is  extremely  close 
to  that  achievable  using  gj.  That  is,  r\{gi.f)  *  r) •(/). 
Of  course,  there  are  situations  where  even  simpler 
detector  nonlinearities  are  desirable.  In  the  next  sec¬ 
tion  we  examine  the  performance  of  three  commonly- 
used  very  simple  suboptimum  nonlinearities. 

ID  THREE  VERY  SIMPLE  SUBOPTIMUM  NON  LINEARI¬ 
TIES 

In  Figure  1  the  locally  optimum  nonlinearity  for  a 
Class  A  noise  PDF  is  plotted  for  itO  (since  /  is  sym¬ 
metric.  gj  is  antisymmetric,  i.e.  gj(-z)  a  -gj(z)). 
The  parameters  (A=0.35,  I"  =  0.0005)  of  the  Class  A 
model  used  in  Figure  1  are  used  In  [9-12]  to  fit 
"interference  (probably)  from  nearby  powerline,  pro¬ 
duced  by  some  kind  of  equipment  fed  by  line"  [9],  The 
vertical  dashed  lines  at  x  *  0.06  and  x  a  0,10  divide 
the  x -axis  into  three  regions  Sj  =  [ lx  |«0.06j, 
St  =  [0.06<|x  i*0.  lOj,  and  5s  *  J  |x  j  >0. 10]  which  are 
the  regions  where  e/(x)  is  approximately  linear  (5|). 
returning  to  zero  (St).  and  approximately  sero  (5g). 
Evaluating  the  probability  under  /  of  each  of  these 
regions  (or,  more  intuitively  the  fraction  of  the  data  we 


h(z) 


h(z) 


h(z)dz. 
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should  expect  to  fall  In  each  region)  we  have  that 
Pr(Sl)  m  0.71.  Pr(St)  *  0.01.  and  Pr(5f)  »  0.28.  Thus 
we  see  that  all  but  about  IX  of  the  time  the  data  will 
fall  in  the  approximately  linear  region  or  the  approxi¬ 
mately  zero  region.  This  leads  us  to  believe  that  gj 
can  be  closely  approximated  by  a  blanker  g%  (also 
called  a  hole  punoher)  which  is  shown  in  Figure  2a.  For 
comparison  we  have  also  examined  the  performance  of 
a  soft  limiter  gh  (or  clipper)  and  a  hard  limiter  gfa  (or 
sign  detector]  which  are  shown  in  Figures  2b  and  2c, 
respectively. 

In  Table  1  we  have  given  the  processing  gain 
achievable  using  the  locally  optimum  nonlinearity  gj, 
the  blanker  gj,  the  soft  limiter  gjL.  and  the  hard  lim¬ 
iter  gm  (Note  that  the  stars  on  ga  and  g 31  indicate 
that  the  optimal  value  of  c  is  used).  We  have  included 
each  of  the  examples  used  in  [l]  as  well  as  two  others. 
In  each  case  we  see  that  the  blanker  is  nearly  optimal 
while  the  soft  limiter  and  the  hard  limiter  have  sub¬ 
stantially  less  than  optimal  performance.  The  one 
exception  to  this  is  the  last  example  (A  =  1.0.  r*0.1) 
where  the  soft  limiter  is  nearly  optimal  and  even  the 
hard  limiter  outperforms  the  blanker.  Not  surprisingly 
the  locally  optimum  nonlinearity  for  this  case  (shown 
in  Figure  3)  is  more  closely  approximated  by  a  soft 
limiter  than  a  blanker.  We  must  stress  though  that, 
based  on  our  experience,  this  seems  to  be  an  unusual 
case.  In  fact  Table  1  is  quite  representative  of  our 
findings  in  general. 

Another  issue  of  importance  when  considering 
various  detector  nonlinearities  is  that  of  robustness  or 
sensitivity.  Since  the  blanker  and  soft  limiter  each 
only  depend  on  one  parameter  (the  "cut-oil"  parame¬ 
ter  c).  it  is  fairly  straightforward  to  examine  their 
robustness.  (Note  that  a  hard  limiter  does  not  depend 
on  such  a  parameter.)  In  Figure  4  we  plotted  the  pro¬ 
cessing  gain  achievable  using  g%  and  g$L  versus  the 
cut-off  parameter  c  for  the  example 
(A  =0.35  r =0.0005)  considered  in  Figure  1.  It  would 
seem  from  the  smoothness  and  flatness  of  the  curves 
in  Figure  4  near  their  respective  maxima  that  both  the 
soft  limiter  and  the  blanker  are  quite  insensitive  to 
variations  in  the  cut-off  parameter. 

On  the  other  band,  in  Figure  5,  the  same  two 
curves  are  plotted  using  a  different  scale  on  the 
abscissa.  This  new  scale  is  not  the  cut-off  parameter  c 
but  the  probability  of  the  set  [-cscxcc]  under  the 
Class  A  PDF  /  (x ).  As  mentioned  above,  this  can  be 
thought  of  as  the  fraction  of  the  data  which  we  can 
expect  to  fall  in  the  linear  region  of  the  detector  non¬ 
linearity  (cf.  5 1  in  the  first  paragraph  of  this  section). 
Since  any  estimate  of  c*.  the  optimal  value  of  the  cut¬ 
off  parameter,  would  presumably  come  from  some  ver¬ 
sion  of  an  empirical  PDF  (see  [12]).  Figure  5  is  likely  to 
be  a  more  reasonable  way  than  Figure  4  to  examine 
the  sensitivity  of  the  blanker  and  soft  limiter  to  uncer¬ 
tainties  in  estimating  c*.  Results  similar  to  Figure  5 
have  also  been  obtained  for  the  Middleton  Class  B 
(broadband)  noise  model  by  Ingram  and  Houle  [IB]. 

It  is  not  unreasonable  at  first  thought  to  assume 
that  this  change  in  scale  would  cause  little  change  in 
the  relative  smoothness  and  flatness  of  the  two  curves. 
In  fact  Figure  5  shows  quite  strlkir^ly  that  the  blanker 
is  very  sensitive  while  the  soft  limiter  is  very  insensi¬ 
tive  near  their  respective  maxima.  This  example  is 
again  quite  representative  of  a  wide  variety  of  other 
cases. 
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nr.  A  ROBUST  APPROXIMATION  SCHEME 

In  tbit  section  we  present  •  new  scheme  for  detec* 
tion  in  noise  which  is  reasonably  well-modeled  by  the 
Class  A  density  (2).  This  scheme  requires  some  back¬ 
ground  on  the  theory  of  asymptotically  robust  detec¬ 
tion.  This  theory  was  first  developed  in  a  paper  by 
Martin  and  Schwartz  [IS]  for  the  case  of  nearly  Gaus¬ 
sian  noise.  It  was  then  extended  by  Kassam  and  Tho¬ 
mas  [20]  to  mixture  models  with  general  "nominal" 
noise  densities.  We  now  present  the  needed  results 
from  [20]. 

In  [20]  an  c-mixture  (or  c-contaminated)  clast  of 
PDFs  of  the  form 

F*  }/  I  /(*)  ■  ( 1  -e)A0(x ) ♦  e A j(x ),  hicHj  (10) 

was  considered,  where  h0  is  a  fixed  element  of  H  satis¬ 
fying  certain  regularity  conditions  and  H  is  the  class  of 
all  bounded  symmetric  PDF's.  They  showed  that  the 
detector  nonlinearity  which  has  the  best  lower  bound 
over  the  class  F  on  asymptotic  performance,  i.e..  the 
minimax  locally  optimum  detector  nonlinearity 
has  the  form 

-ho(-c) 

— r- 7 - r—  *<-e 

fco(-c) 

<<'’■  i*i“  <“> 


where  c  is  chosen  to  satisfy 

/*•<*>“ <«> 

As  we  see  from  (11).  the  robust  detector  non¬ 
linearity  is  equal  to  the  locally  optimal  nonlinearity  for 
the  nominal  PDF  hB  (see  (8))  for  |z  I  less  than  the  cen¬ 
soring  constant  c,  and  it  is  constant  for  |x|  greater 
than  c.  i.e..  is  a  censored  version  of  g^ 

In  [1]  we  saw  that  in  most  cases  the  locally 
optimum  detector  nonlinearity  g}  for  the  Class  A  PDF 
(2)  can  be  closely  approximated  by  the  locally 
optimum  nonlinearity  gi  for  the  PDF  containing  just 
the  first  two  terms  (i.e.,  /•**  as  in  (5)).  In  the  cases 
where  there  was  any  substantial  difference  between 
the  two  nonlinearities  it  was  in  the  tails.  Beyond  4  or  5 
standard  deviations,  the  two  term  nonlinearity  g{  con¬ 
tinued  its  linear  path  while  the  infinite  term  nonlinear¬ 
ity  flattened  virtually  to  the  horizontal.  In  Figure  6 
these  two  nonlinearities  have  been  plotted. 

The  horizontal  (dotted)  line  added  to  Figure  8 
corresponds  to  the  robust  detector  nonlinearity  for 

the  class  F  with  Ao  =  /(,)  >nd  t  =  J]  That  is. 

instead  of  simply  discarding  all  the  terms  in  the  Class 
A  PDF  beyond  the  first  two  as  in  [l],  we  have 
"robustified"  with  respect  to  this  remainder.  It  is  clear 
from  Figure  6  that  this  robustified  two  term  nonlinear¬ 
ity  g!  may  often  be  an  even  better  approximation  to 
g}  than  gi  Our  conviction  has  been  that,  in  those 
cases  where  the  two  term  locally  optimum  nonlinearity 
gi  does  not  perform  in  a  nearly  optimal  fashion 
(assuming  the  noise  actually  has  the  Class  A  density  / 
as  in  (2)).  the  robust  version  g{  will.  However,  as  of 
this  writing  we  have  been  unable  to  find  a  single  case 
where  gi  does  not  perform  virtually  as  well  as  gj. 


Still  there  are  justifications  for  this  robust 
scheme.  First,  in  every  case  we  have  considered  the 
robuet  version  performs  about  as  well  or  even  better 
than  the  two  term  optimal  version.  Even  though  these 
differences  are  negligible  they  do  not  contradict  the 
conviction  stated  above.  Of  course,  this  requires 
further  exploration.  Second,  it  has  been  shown  by 
Berry  [13]  that  the  Class  A  model  would  generaly  be 
useful  when  the  noise  is  dominated  by  two  com¬ 
ponents:  receiver  noise  and  one  strong  interferer.  If 
this  is  the  case,  the  Class  A  model  may  be  chosen 
solely  on  the  basis  of  its  fit  for  the  first  two  (dominat¬ 
ing)  terms  and  the  remaining  terms  may  not  really  be 
accurate.  Clearly,  due  to  its  minimax  design  the 
robust  version  will  be  lees  sensitive  to  these  inaccura¬ 
cies  than  the  two  term  optimal  detector.  In  fact,  it 
might  even  outperform  the  "optimal"  nonlinearity  0/  if 
the  true  npise  has  other  small  interferers  which  are 
not  at  the  appropriate  distances  and  powers  necessary 
to  be  accurately  described  by  the  Class  A  model 
beyond  the  first  two  terms. 

V.  SUMMARY 

In  this  paper  we  have  examined  the  performance 
of  some  simple  suboptimum  detector  nonlinearities. 
For  most  Class  A  noises  the  blanker  has  nearly  optimal 
performance  while  the  soft  limiter  and  bard  limiter 
have  significantly  lower  performance.  On  the  other 
hand  the  performance  of  the  blanker  seems  to  be  far 
more  sensitive  to  errors  in  estimating  the  optimal 
cut-off  parameter. 

In  [l]  we  showed  that  the  Middleton  Class  A  noise 
model  can  often  be  approximated  closely  by  the  (nor¬ 
malized)  sum  of  just  the  first  few  terms.  In  fact,  in 
many  cases,  two  terms  are  sufficient.  This  was  espe¬ 
cially  clear  when  we  looked  at  the  efficacy  of  the 
detector  nonlinearity  which  is  locally  optimum  for  two 
terms  of  the  Class  A  model  and  found  it  comparable  to 
the  efficacy  of  the  full  locally  optimum  detector.  In 
this  paper  we  also  developed  a  robust  scheme  which 
needs  some  further  examination  but  which  sre  feel 
offers  a  number  of  advantages  over  the  methods  of  [1]. 
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Figure  2a.  The  blanker  (hole  puncher). 


Figure  2b.  The  soft  limiter  (clipper). 


Figure  2c.  The  hard  limiter  (sign  detector). 


Table  1.  Processing  gain  (efficacy)  achievable  using 
optimal  (f  *).  blanker  \g}).  soft  limiter  (p«.).  and  hard 
limiter  (p*,)  nonlinearities 


Figure  1.  Locally  optimum  detector  nonlinearity  for 
Class  A  noise  PDF  /  (linear  scales)  (A *0.36,  P*0  0006) 
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n*w*  9.  Locally  optimum  detector  nonlinearities  gi  Figure  5.  Efficacy  (processing  gain)  of  soft  limiter  pj 

and  if  and  robust  nonlinearity  g  f  (dashed  line)  and  blanker  p|  (<4=0  36,  r=0.0005).  Abscissa  seal 

(A  =0.36,  r =00005)  different  from  Fig.  4. 
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